Collisional properties of sympathetically cooled 39 K 
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We report the experimental evidence of the sympathetic cooling of 39 K with 87 Rb down to 1 /iK, 
obtained in a novel tight confining magnetic trap. This allowed us to perform the first direct mea- 
surement of the elastic cross section of 39 K below 50 /iK. The result obtained for the triplet scattering 
length, <2t = —51(7) Bohr radii, agrees with previous results derived from photoassociation spectra 
and from Feshbach spectroscopy of 40 K. 
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I. INTRODUCTION 

In the field of ultracold and quantum gases Potas- 
sium gained a key role since it is the only alkali, be- 
side Lithium, for which a stable fermionic isotope exists. 
Furthermore, since the sympathetic cooling of fermionic 
40 K with 87 Rb, the most widespread atom in the labora- 
tories devoted to laser cooling, is particularly favorable, 
it is not surprising that so many experiments with ul- 
tracold fermions were carried out or are planned with 
this atom. In contrast very little attention was paid to 
the two other (bosonic) isotopes of Potassium, 39 K and 
41 K. These isotopes offer the possibility of creating dou- 
ble species Bose-Einstein condensates, which display rich 
and interesting quantum phase diagrams when trapped 
in optical lattices p], 0, [3, 0]. After the seminal work 
where a double species 41 K- 87 Rb BEC was produced in 
a magnetic trap [5[, essentially no experiments have ad- 
dressed these topics. 

Along these guidelines, we started an experiment de- 
voted to the exploration of degenerate Bose-Bose mix- 
tures in optical lattices. In this work we report that 
also 39 K can be sympathetically cooled with 87 Rb al- 
though this is much more difficult due to its unfavorable 
laser cooling and scattering properties, thus demonstrat- 
ing that all the stable isotopes of Potassium can be cooled 
to ultralow temperatures. 

Laser cooling of 39 K has several disadvantages if com- 
pared to alkalis such as 87 Rb or even 40 K. These dis- 
advantages are related to the hyperfine structure of the 
P3/2 level whose separation is of the same order of mag- 
nitude as the natural linewidth of the cooling transition 
\F = 2,m F = +2) \F = 3,m F = +3). The result- 
ing strong optical pumping towards the \F = 1) manifold 
must be countered by a repumping light with an intensity 
comparable to that of the light driving the cooling transi- 
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tion. A careful balance of the intensity and the detuning 
of these two laser frequencies, unsuitable for optimized 
loading of a Magneto- Optical- Trap (MOT), is the only 
way to obtain temperatures at the Doppler limit by laser 
cooling [6[. Furthermore, the s-wave elastic cross section 
of 39 K at zero temperature is more than one order of mag- 
nitude smaller than that of 87 Rb and, due to the attrac- 
tive character of the interaction, the Ramsauer-Townsend 
minimum [7] occurs at a temperature of around 320 /iK 
where the contribution of other partial waves is still small 
[8[. All these features strongly hamper the efficiency of 
evaporative cooling. Note that, due to Pauli blocking, an 
even more dramatic decrease in the evaporative cooling 
efficiency occurs for the fermionic 40 K. 

For these reasons the most widely employed experi- 
mental technique to reach ultralow temperatures in a K 
gas is sympathetic cooling with 87 Rb. This amounts to 
producing a mixed MOT of the two species and then 
evaporatively cool 87 Rb allowing K to thermalize with it. 
With this technique it has been possible to reach quan- 
tum degeneracy of 40 K @ and 41 K [10], but no attempts 
have been made to test this method on 39 K. For this 
mixture in fact, the inter-species cross section is again 
more than one order of magnitude smaller than the one 
of 87 Rb- 40 K and 87 Rb- 41 K [llDul. 

We demonstrate here for the first time that 39 K can be 
sympathetically cooled to 1 /iK. Therefore the work pre- 
sented in this paper is a critical step towards the produc- 
tion of a Bose-Einstein Condensate (BEC) of 39 K. Such 
a system is a promising candidate for the realization of a 
BEC with interaction tunable around zero, since a broad 
Feshbach resonance is predicted around 400 G [HEl. 

The achievement of sympathetic cooling allowed us to 
perform the first direct measurement of the s-wave elastic 
cross section of 39 K. Assuming from earlier work on pho- 
toassociation spectra [13] and on Feshbach spectroscopy 
of 40 K [14[ the attractive character of the interatomic 
interaction (i.e. the negative sign of the s-wave scatter- 
ing length), our results agrees with [H, Q within 1.7 
combined standard deviations. 

The time required for sympathetic cooling depends 
crucially on the collision rate between the two species. 
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If this time is of the same order of magnitude as the 
lifetime of the sample, the density will drop during the 
evaporation and cooling will eventually stop. Since the 
cross section for 39 K- 87 Rb collision is smaller than in the 
case of 40 K and 41 K, one has to increase the lifetime of 
the sample and/or the confinement before starting sym- 
pathetic cooling. In the experiments reported here, the 
confinement of the 39 K sample is increased by more than 
a factor 2 with respect to [5[ employing a new kind of 
magneto-static trap of size intermediate between micro- 
traps, where microscopic current-carrying wires are laid 
on a chip, and traditional magnetic traps created with 
multiple winding coils [151 ]. 

The outline of the paper is the following: in section 
[II] we will present a brief description of the experimental 
apparatus and the results about sympathetic cooling of 
87 Rb and 39 K. In section!]]] we will introduce the exper- 
imental procedure for studying ultracold collision in 39 K 
and we will present the theoretical analysis of the data. 
Finally, we will conclude with the perspectives of 87 Rb- 
39 K evaporative cooling and of a 39 K BEC. 



II. EXPERIMENTAL APPARATUS 
A. Laser system and MOT loading 

In order to provide the maximum flexibility for the 
loading of the mixed MOT, our experimental apparatus 
features two distinct two-dimensional MOT's (2D-MOT) 
for the separate precooling of 87 Rb and 39 K. A complete 
description of our laser system and a thorough charac- 
terization of our 39 K 2D-MOT can be found in Ref. jig . 
The 2D-MOT employed for 87 Rb has exactly the same 
structure and similar performances. In the experiments 
we used 130(110) mW for the cooling light and 40(6) 
mW of repumping light for 39 K ( 87 Rb). The atomic flux 
provided by the two 2D-MOT's is collected inside a mixed 
MOT formed in a ultra high vacuum chamber. The MOT 
is formed by a standard configuration with six indepen- 
dent laser beams with 25 mm diameter and a total power 
of 45 (40) mW for the cooling light and 24 (6) mW for the 
repumping light of 39 K ( 87 Rb). These beams are split 
from the output of one single mode fiber in which we 
inject the four different frequencies needed for the two 
species, thereby guaranteeing a perfect geometric over- 
lap of the cooling and repumping light and of the two 
MOT's. The magnetic field gradient of about 16G/cm 
used in the MOT is provided by a pair of anti-Helmholtz 
coils operating at 4 A. 

We typically load around 2 x 10 9 atoms of 87 Rb and 
5 x 10 6 atoms of 39 K in 10 and 1 s respectively. We can 
control very precisely the number of 39 K atoms loaded 
in the MOT adjusting the time of operation of the 2D- 
MOT. After the loading of the MOT we reduce the re- 
pumping intensity and increase the magnetic field gra- 
dient to compress the cloud, we apply a 3.5 ms stage of 
polarization-gradient cooling, optically pump the atoms 



of both species towards the \F = 2, mp = +2) state and 
then trap them in a purely magnetic trap. 



B. Magnetic trap and evaporation 

Our magnetic trap represents a compromise between 
the usual magnetic traps generated by coils placed out- 
side the vacuum system, typically low confining (^ max ~ 
100 Hz), high power consuming (~ lkW), several cen- 
timeters in size, which produce the largest BEC's (> 10 6 
atoms), and the so called microtraps operating in vac- 
uum, providing a high confinement (^ ma x ~ 1 kHz), re- 
quiring low electrical power (~ 1W), that easily fit on 
a microchip, but usually produce smaller condensates 
(~ 10 4 atoms). 

The potential generated by our trap is of the Ioffe- 
Pritchard type and the structure has both current con- 
ductors laying on a chip and free standing [15[ . The free 
standing conductors are built from a single oxygen-free 
copper tube machined to form the four Ioffe bars and two 
partial rings that provide the end caps of the potential in 
the axial direction. The axial confinement is increased by 
a circular copper trace laid on a direct-bond copper chip 
that interfaces the current leads and the free standing 
conductors. Electrical and mechanical contact between 
the parts is obtained by vacuum brazing and the current 
leads are brought outside the Ultra-High Vacuum envi- 
ronment by a custom-made electrical feedthrough. To 
allow optical access along the horizontal axial direction, 
the feedthrough is hollow and terminates with a vacuum 
viewport; optical access in the radial direction is easily 
obtained through the 2 mm gaps between the Ioffe bars. 
Since these gaps are too small to allow the efficient pro- 
duction of a MOT, the magnetic trap, hereafter called 
millimetric-trap or mTrap, is displaced by 26 mm from 
the center of the MOT. A big coil placed outside vac- 
uum is used to adjust the bias field of the trap thereby 
changing the radial confinement. The confinement pro- 
vided by our mTrap at 95 A current and 6 G bias field is 
uo r = 2tt x 447(7) Hz and u z = 2tt x 29.20(1)Hz for 39 K. 

The displacement of the atoms from the MOT posi- 
tion to the mTrap is accomplished holding them inside a 
quadrupole trap which is moved with a motorized trans- 
lation stage. The quadrupole trap is formed ramping the 
current of the MOT coils up to 65 A and the motion 
takes less than 500 ms in order to reduce the losses due 
to Majorana spin-flip. 

While the magnetic field of the mTrap adds to that of 
the transfer trap along the radial directions, along the 
axial direction the mTrap generates a field directed in 
the opposite direction so that an adiabatic transfer of 
atoms would be only possible decreasing the trap depth 
during the process and therefore reducing the number 
of transferred atoms. To overcome this potential pitfall 
we employ a mixed approach to load the trap: we ramp 
up the current of the mTrap in 150 ms, a time which is 
fully adiabatic only for the radial direction. The overall 
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III. CROSS-THERM ALIZATION 
MEASUREMENT 



A. Introduction 



FIG. 1: (Color online) Experimental sequence showing sym- 
pathetic cooling between 87 Rb (left) and 39 K (right). All the 
images are taken after 1 ms expansion. The phase space den- 
sity at 1 /xK is 0.02 and 0.01. Images are taken along the 
vertical radial direction. 



efficiency of the transfer from the MOT to the mTrap 
is around 30%, and we typically trap 7 x 10 8 87 Rb and 
4 x 10 5 39 K atoms. 

Once atoms are loaded inside the magnetic trap we per- 
form evaporation on the 87 Rb cloud by transferring the 
hottest atoms to the antitrapped \F = l,mf = +1) state 
with a microwave sweep around 6.8 GHz. In addition we 
apply a second ramp to remove unwanted atoms from the 
\F = 2,mf = +1) level. The 39 K is unaffected by these 
sweeps, owing to the huge difference in frequency with 
the corresponding hyperfine transition (460 MHz), and 
only minor losses occur during the evaporation. Provided 
that the number of 87 Rb atoms Nr^ at a given temper- 
ature is still larger than the initial number of 39 K atoms 
Nk, the former can act as coolant establishing thermal 
equilibrium, as it is shown in Fig. [TJ For this reason the 
initial number of 39 K atoms is reduced in the pictures 
from 2.7 x 10 5 in the first picture to 5.9 x 10 4 of the last 
one at a temperature of 1.0(1) fiK. 

We experimentally found that a sufficient condition to 
reach thermal equilibrium is Nr^/Nk > 2.7, but no at- 
tempt has been made to measure the inter-species elastic 
cross section since it can be obtained with high precision 
from the results of Feshbach spectroscopy of the 40 K- 
87 Rb mixture We used instead the known value 

dK-Rb = +36ao to guide the optimization of the evap- 
oration. Indeed we experimentally verified that in pres- 
ence of 39 K atoms, in order to optimize the number of 
87 Rb atoms below a few /iK, the evaporation ramp must 
be slowed down in the last part to be sure that the two 
species are always in thermal equilibrium. 



We measure the cross section of 39 K by analyzing the 
relaxation of the atomic cloud towards thermal equilib- 
rium. The technique can be described as follows. The 
dynamics of an ideal gas isolated and confined in an 
harmonic potential is completely separable meaning that 
the different spatial degrees of freedom are decoupled. 
If such a sample is prepared with different energy dis- 
tributions along different directions this difference will 
not vary with time. On the other hand, due to inter- 
actions, the collisions drive the system towards thermal 
equilibrium, namely a state in which the total energy is 
conserved and the effective temperatures of the differ- 
ent degrees of freedom are equal. If the initial state has 
different effective temperatures, the relaxation towards 
thermal equilibrium can provide direct information on 
the collisions. In our system this is accomplished in the 
following way: after obtaining a cloud of cold 39 K as de- 
scribed in Section HI 8T Rb is blowed away from the trap 
with a pulse of resonant light and the radial degrees of 
freedom of K are excited by parametric excitation, ob- 
tained modulating the value of the trap bias field at the 
frequency of radial confinement for 100 ms. Due to our 
elongated trap geometry this excitation is highly selec- 
tive and the axial degree of freedom maintains its effec- 
tive temperature. Following this excitation we allow the 
cloud to relax for a time t w after which we image the 
cloud along the vertical radial direction by means of in 
situ absorption imaging. From the measured widths w 
of the Gaussian density profile we can calculate the av- 
erage potential energy along the radial (r) and axial (z) 
direction as 

Ei = ^k B Ti = ^(uiWi) 2 i = r,z . (1) 

where m is the atomic mass. Repeating this measure- 
ment for different t w allows us to measure the relaxation 
of the ratio between the axial and radial temperatures. 
As shown in Fig. [2] this is well fitted by an exponential 
decay. The time constant r of this decay is related to the 
elastic scattering cross section and to the density of the 
sample by the following relation: 

r _i = Tel = n(av) 
a a ' 

where a is the number of collisions required to an atom to 
reach thermal equilibrium, j e i is the intraspecies collision 
rate for 39 K, n is the average density, a is the scattering 
cross section, v is the relative velocity of two colliding 
atoms and (•) indicates averaging over the velocity dis- 
tribution in the cloud. 

If one indicates the geometric average of the tempera- 
ture along the different direction as T = (T^T^) 1 / 3 and 
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introduces an average, temperature-dependent cross sec- 
tion defined as: 



a(f) = 



(v) 



since the average density of a harmonic trap is: 



US?) Nk ^ 



and the average 
4 (for) 1 ' 2 ( 7rm )-i/2 > 
pressed as: 



relative velocity is (v) 
the relaxation rate can be 
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We will come back to the relation between a and the 
actual energy-dependent cross-section in Sec. IIII CI 

Eq. ([5j) holds only if one assumes that the average den- 
sity of the sample does not change during the relaxation 
process. Actually the loss rate induced by collisions with 
atoms of the background gas is about 0.03 s _1 for the ex- 
periments reported in this work. If this rate is not negligi- 
ble as compared to the measured relaxation rates, it must 
be taken into account: one expects that the reduction of 
the density slows down the relaxation so that after an 
initial decay with a rate corresponding to Eq. ([5]), equi- 
librium apparently sets for a value of E z /E r lower than 
1. This picture is further complicated by the unavoidable 
heating rate present in the experiment that increases the 
temperature throughout the measurements: in our sys- 
tem we measure a different heating rate in the axial and 
radial direction so that it is possible to achieve equilib- 
rium also for a value of E z /E r greater than 1. 

In principle, one could take these effects into account 
having the equilibrium value of the energy ratio as a free 
parameter: for our measurements however the equilib- 
rium value is almost always consistent with 1. We there- 
fore choose to fix it to 1 in order to improve the esti- 
mation of the time constant, as it is shown in Fig. [2j 



B. Measurements of the elastic cross section 

In order to reduce systematic errors on the determi- 
nation of <j, we performed several measurements at dif- 
ferent densities and for two different temperatures of 
29 and 16 jaK. Furthermore, to check that the inhomo- 
geneous magnetic field present in our in situ imaging 
was not a source of error, the dataset at 16 fiK is taken 
with a slightly different procedure. First, we adiabati- 
cally decompress the trap to uo r = 2tt x 290(1) Hz and 
uj z = 2tt x 21.24(1) Hz and then blow away 87 Rb, apply 
parametric heating at cj r , wait for a variable time and 
image the cloud after a 2 ms expansion. 

The results of our measurements are shown in Fig. [3] 
where we report the relaxation rate 7 = r _1 as a function 
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FIG. 2: (Color online) Plot of the relaxation dynamics of 
an ultracold sample of 39 K after parametric heating in the 
radial direction. Data are taken after expansion (see text). 
Each point is an average of several experimental realization 
and the line is an exponential decay fit with equilibrium value 
fixed to 1, initial value and time constant as free parameters. 
Average initial number of atoms is 379(10) x 10 3 and initial 
average temperature is 16.2(7) /nK. From the fit we obtain 
r = 3.66(46) s. 



of the initial number of atoms Nk- Figure [31(a) shows the 
measurements taken after expansion with an average ini- 
tial temperature T = 16(1) /iK, while Fig. [31(b) presents 
measurements taken in situ with T = 28.9(1.3) /jK. One 
can see that the relaxation rate has the expected linear 
dependence on the number of atoms and the extrapola- 
tion towards zero is consistent with zero. 

In a Ioffe-Pritchard trap, a sufficiently big atomic cloud 
can experience a potential that is not strictly separable 
so that different degrees of freedom are coupled and re- 
laxation can occur even in the absence of collisions. We 
refer to this process as ergodic mixing. A priori, for equal 
harmonic frequencies, ergodic mixing could play a more 
prominent role in our mTrap than in the usual traps due 
to its small size. For this reason we take ergodic mixing 
into account by separating the component of the relax- 
ation rate, which is linear in TV, from the extrapolation in 
the limit of zero density where relaxation can only occur 
through ergodic mixing. One can therefore assume that 
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where imix is the time constant of the ergodic mixing 
process. The two parameters A and jrnix are obtained 
from a linear fit of the experimental data, as shown in 
Fig. [3l By using Eq. ([5]) one has 
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As pointed out above, since the values for jmix are con- 
sistent with zero in both the datasets, we can conclude 
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FIG. 3: (Color online) Plot of the measured relaxation rate 
as a function of the number of atoms in the sample after ex- 
pansion (a) and in situ (b). The measured initial average 
temperature is 16(1) /xK and 28.9(1.3) /iK respectively. The 
slope of the linear fit to the data is related to the elastic 
cross-section, while the intercept is a measure of ergodic mix- 
ing as described in the text. The fit results are: jmix = 
-0.003(0.065) s" 1 and A = 5.5(1.6) x 10" 7 s _1 for figure (a) 
and <y m ix = -0.0005(0.034) s" 1 and A = 4.0(1.0) x lO^s" 1 
for figure (b). Error bars on the vertical direction are statis- 
tical error on the exponential decay fit, while on the horizon- 
tal direction they are the statistical fluctuation of the initial 
number of atoms not including the calibration uncertainty. 



that ergodic mixing plays a negligible role in our mea- 
surements. 

The resulting values for the temperature-dependent 
average elastic cross section, computed from Eq. ([7]) 
with a = 2.7, as obtained from numerical simulation 
described in Sec. UTTPl are 2.2(0.8) x l(T 12 cm 2 and 



0.91(0.22) x 10- 12 cm 2 for 16 and 29 /iK respectively. 
These values are a very good approximation of the ac- 
tual cross section as we show in the next section. 



C. Extracting the value of the scattering length 

As pointed out above, the measured cross sections de- 
pend on the temperature. In order to extract informa- 



tion on the interatomic potential, one has to make some 
assumption on the behavior of the cross section as a func- 
tion of temperature. A simple partial wave expansion of 
the scattering amplitude shows that, for the experimen- 
tal condition reached in this work, collisions can only 
occur at I = angular momentum, odd values being sup- 
pressed on symmetry ground and even values due to the 
low temperature [27[. In a very general form the s-wave 
cross section for identical boson is: 



a(k) = 87T|/ (fc)| s 



(8) 



where k = mv/(2h) is the relative wavevector and fo(k) 
is the s-wave scattering amplitude that can be expressed 
in terms of the S-matrix phase shift 5o(k): fo(k) = 
e iS °( k ) S m5 (k)/k. 

In the limit of vanishing energy, given that the poten- 
tial decays more rapidly than 1/r 3 , we can express the 
cross section as a function of a single parameter, the s- 
wave scattering length which, for two alkali atom in the 
stretched \F = 2,mi? = +2) state, is the triplet one ar ' 



\fo(k)\ 2 
a(k) 



(9) 
(10) 



For higher temperature one has to calculate the scat- 
tering amplitude up to order k 2 : fo(k) — (— + 
ik + \k 2 r e + . . .) _1 , where r e denotes the effective range, 
which is a function of the scattering length ot and the 
Cq coefficient of the van der Waals potential [17| 



a(k) 8tt 



(1 - \k 2 r e a T ) 2 



k 2 cl 2 ^ 



(ii) 



Since the Cq coefficients are well known for all alkali 
dimers [l8|, Qjl, we numerically invert Eq. (j3]) after in- 
serting Eq. (fTTj), averaging over the Boltzmann thermal 
distribution and assuming ar < 0, to obtain the scat- 
tering length: ar = — 67(ll)ao from data at 16 /iK and 
a T = — 48(5)ao from data at 29 /iK [28]. For comparison, 
we notice that, if we neglect the effective range, these 
values would be — 57(11) ao and — 36(5)ao, respectively. 
In the first case the two values differ by approximately 
1.57 combined standard deviations, while in the second 
case the difference is 1.75. The quoted uncertainties de- 
rive mainly from the fits of Fig. [3] and the atom number 
calibration (±20%), which was done independently for 
each data set. Finally, we take a weighted average of 
the two scattering length values calculated with effective 
range and multiply the associated uncertainty by a factor 
\[\ 2 = 1.65, to set the confidence level to 68% [20]: the 
result is ar = — 51(7)ao- 

Prior to this work, the 39 K triplet s-wave scatter- 
ing length has been measured by photoassociation spec- 
troscopy, <2t = — 33(5)ao [HI, and by mass-scaling from 
rethermalization of 40 K, ar = — 37(6)ao [14J- Indeed our 
measurement is the first direct determination of the elas- 
tic cross section for this isotope. To different extents, all 
reported values of the scattering length depend on the 
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Cq coefficient of the van der Waals long range potential: 
we take C 6 = 3927 a.u. while in Ref. 0,13 C 6 is 
assumed equal to 3897 a.u. [18]. We checked that our 
value changes with Cq with a rate S\a t \ = O.OOlao x SCq 
[a.u.], approximatively 50 times smaller than in Ref. [13[ 
(no such rate is available in Ref. [lj]) and therefore the 
difference in the Cq does not significantly change the de- 
rived scattering length. As our result differs from the av- 
erage of the published values of 1.70 combined standard 
deviation, we conclude that the agreement is satisfactory. 



D. Numerical simulations 

As shown in Eq. Q, to compute the value of the scat- 
tering length from the fit to the data one needs to know 
the value of the parameter a, which is the average number 
of collisions per particle needed for thermalization. This 
parameter depends on the temperature and the trap fre- 
quencies in a non trivial way and can vary from 2.4 to 3.4 

[Hi. Following mmim we estimate this parameter 

from a numerical simulation of the system. Taking ad- 
vantage of the low density and small atom number of the 
sample, we make a direct simulation of the gas in which 
we consider 3D position and velocity of every single atom. 
Choosing a time step St much smaller than both the aver- 
age collision time 7^ and the faster timescale of single 
particle dynamics, namely a;" 1 , one may assume that 
during the time interval St interactions between atoms 
are weak enough that they decouple with the center of 
mass motion and that the force experienced by the atoms 
varies little. Therefore the dynamics of the gas can be 
directly simulated with the following procedure: 

1. Position and velocity are updated according to ex- 
ternal forces with a Verlet integrator. 

2. Every M steps, the positions of the atoms in real 
space are discretized on a lattice of spacing Sx and 
if two atoms are located at the same site a collision 
test is made. 

3. If the collision test is positive, the collision is re- 
solved using simple classical mechanics and s-wave 
scattering: 

! IvJ ~ 

v 2 vcm - ^ e R 

where indicates the velocity of particle i after the 
collision, vcm and v r are the center of mass and 
relative velocity respectively, before the collision, 
and e# is a random direction on the unit sphere. 

The collision test consists in comparing a random real 
number y, uniformly distributed in [0, 1), with the colli- 
sion probability calculated according to kinetic theory: 

p = a\v r \MStSx~ 3 (12) 
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FIG. 4: (Color online) Plot of the simulated energy ratio 
(black dots, left scale) and collisions per particle (blue tri- 
angles, right scale) as a function of time for 20 x 10 3 atoms 
in a harmonic trap whose frequencies are 440 and 29.4 Hz in 
the radial and axial direction respectively. The red curve is 
an exponential decay fit. The value a = 2.72 can be obtained 
directly from the graph identifying the time constant of the 
decay (7.74 s). 



where a is the elastic scattering cross section and \v r \ is 
the modulus of the relative velocity; collision is processed 
if y < p. In order to reproduce a correct scattering rate, 
one must choose Sx so that the average number of atoms 
in a volume Sx 3 is much smaller than 1 and M so that 
the distance traveled by an atom during the time M St is 
of the order of Sx. 

In Fig. [4] we report a typical result of the simulation 
showing the radial to axial energy ratio and the number 
of collisions per atom as a function of time. The simula- 
tion is performed in a cylindrical harmonic trap with ra- 
dial and axial frequency of 440 and 29.4 Hz respectively. 
The number of atoms is 2 x 10 4 and the initial energy 
distribution has a width of 15(19.5) /uK along the axial 
(radial) direction. The cross section is 8.6 x 10 _16 m 2 so 
that the scattering rate in the simulation is close to the 
experimental one. The value of a = 2.7 can be extracted 
directly from the graph as explained in the caption of 
Fig.H 



IV. CONCLUSIONS 

In summary we have demonstrated sympathetic cool- 
ing of 39 K with 87 Rb, in spite of the low interspecies 
cross section, and we have measured the cross section for 
elastic collision between two 39 K atoms in a pure triplet 
state. We have hence obtained the triplet s-wave scatter- 
ing length <2t , assuming its negative sign and an effective 
range approximation for the s-wave scattering amplitude 
at low but finite temperature. 

The possibility of sympathetic cooling and the knowl- 
edge of the collisional properties of bosonic potassium 
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open the way to its use as a quantum degenerate gas. 
Bosonic potassium isotopes, 39 K and 41 K, are predicted 
to feature Feshbach resonances, several Gauss wide, at 
moderate fields (< 1 kG). As a consequence, optically 
trapped bosonic potassium appears a suitable candidate 
to realize a condensate with a scattering length tunable 
around zero. Such a condensate would be interesting for 
inter ferometric purposes, as the interaction energy limits 
the accuracy of interferometers employing Bose-Einstein 
condensates. In addition, wide Feshbach resonances 
make potassium condensates particularly attractive to re- 
alize double condensates in optical lattices, for quantum 
simulation purposes. For 39 K, in particular, Bohn and 
coworkers have predicted a Feshbach resonance around 
40 G for two atoms in the \F = 1, mp = — 1) state, with 
a width of several tens of Gauss [25[, while another 60 
G broad resonance has been predicted around 400 G for 
atoms in the absolute ground state \F = 1, tuf = 1) [26[. 



While a double condensate of 41 K and 87 Rb was produced 
in a magnetic trap [5], 39 K has not been brought to de- 
generacy so far. This work represents a decisive step in 
this direction as we reached a phase-space density of 0.01, 
with a gain of several orders of magnitude with respect 
to the state-of-the-art 
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